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We study dissipation in a dilute Bose gas induced by the motion of a macroscopic object. A blue- 
detuned laser beam focused on the center of a trapped gas of sodium atoms was scanned both above 
and below the BEC transition temperature. The measurements allow for a comparison between the 
heating rates for the superfluid and normal gas. 

I. INTRODUCTION 

Superfluid flow is a manifestation of quantum mechanics at the macroscopic level. Phenomena like dissipationless 
flow and persistent currents can be traced back to the existence of a macroscopic wavefunction. The gradient of its 
phase gives the superfluid velocity However, superfluid flow is only stable below a critical velocity. Above this 
velocity, excitations of the fluid can be generated, which can be either phonons, rotons (in liquid ''He) or vortices. 

The realization of Bose-Einstein condensation in dilute gases has created a new test ground for the theory 
of quantum fluids. Early experiments confirmed the microscopic foundation of the phenomenon of superfluidity: 
the phonon nature of low-lying collective excitations was observed and the coherence or macroscopic phase 

of the condensate ||^. This picture became more complete during the last year with the observation of vortices 
transverse excitations [0, evidence for a critical velocity and suppression of coUisions [^. Although 
the metastability and therefore limited lifetime of the quantum gas and its mesoscopic size may prevent spectacular 
observations of persistent flows, there is already a clear picture emerging of the many facets of superfluidity in a gas. 

In this article we present measurements which indicate a regime of frictionless flow for a gaseous Bose-Einstein 
condensate. We emulate the motion of a macroscopic object by scanning a focused, blue-detuned laser beam through 
the gas at different velocities. By taking measurements on nearly pure condensates as well as purely thermal ensembles, 
we can compare the relative role of dissipation in both phases. This extends our previous work | pd] , p^ where we studied 
dissipation in a Bose condensed system. 

II. THE EXPERIMENTAL SET-UP 



The experiments were performed in a new apparatus for the production of Bosc-Einstcin condensates of sodium 
atoms. The apparatus has been functioning since January 1999 and incorporates some improvements with respect 
to the original machine of the MIT group. In particular, a Zeeman slower with magnetic field reversal ("spin-flip" 
Zeeman slower) delivered a flux of typically 10^^ atoms s^^. This resulted in trapping 2 — 3 x 10^" atoms in a 
dark-spontaneous force optical trap at a temperature of ~ 1 mK. After 3-4 ms of polarization gradient cooling, atoms 
at a temperature of 50-100 /xK were loaded into an loffe-Pritchard magnetic trap. The trap consisted of four loffe 
bars and two elongated pinch coils, symmetrically located around a quartz cell, combining tight confinement and 
excellent optical access on four sides of the glass cell. Under the tightest confinement, the cylindrically symmetric 
confining potential has radial and axial trapping frequencies Vr = 547 Hz and = 26 Hz. The initial stage of 
evaporative cooling (18 s) proceeded using this trapping configuration; however, the final stage evaporation appeared 
to be limited by inelastic losses at high densities. Therefore, the magnetic confinement was reduced in the last 5-7 
s, finally producing condensates of 2-5x10^ atoms in traps of radial frequencies = 40 — 80 Hz and axial frequency 
Vz = 20 Hz. Since decompression occurred during the evaporation and not after the condensate is produced, as in 
previous work non-adiabatic effects which lead to heating were eliminated. Thus nearly pure condensates (> 90% 
condensate fraction) were produced regardless of the final trapping geometry. 
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FIG. 1. Schematic view of the experimental setup. The condensate was confined in a glass cell using an lofi'e-Pritchard 
magnetic trap. The blue-detuned laser beam was sent through a two-axis AOM scanner and then focused on the condensate 
using the lens LI. The probe light was overlapped using a small mirror and, after passing through the condensate, was focused 
with the lens L2 and sent to the CCD camera. The filter F was used to avoid exposure of the CCD camera to the high intensity 
blue-detuned laser beam. 

A 514 nm laser beam used to stir the atoms was split off from the argon-ion laser that pumped the dye laser 
used for the cooling and trapping beams. Such focused blue-detuned beams have been used previously to eliminate 
Majorana losses in the center of a quadrupole magnetic trap Q and to excite sound waves A single mode optical 
fiber significantly improved the stability with respect to previous experiments ]Tl| . The position of the laser focus 
did not change by more than a few /xm per day and required only minimal adjustment to keep it centered on the 
atomic cloud. The laser was focused to a Gaussian 1/e^ beam diameter of 10 to 13 ^m. The repulsive optical dipole 
force expelled the atoms from the regions of highest laser intensity, and we typically used a ratio of barrier height to 
chemical potential (or temperature, in the case of purely thermal clouds) of 5-10. The condensate had a mean-field 
energy of 60-140 nK, while the thermal cloud temperature was varied between 1 and 6 fiK. Thus, about 200-400 /xW 
of laser power were required for stirring a condensate, and 6-10 mW for stirring a thermal cloud. 

The laser was focused on the center of the cloud by observing the "pierced" condensate using in-situ phase-contrast 
imaging The focus was scanned back and forth along the axial direction using an acousto-optic deflector. The 
velocity of the scan was varied by adjusting the frequency / and keeping the amplitude a fixed, yielding a velocity 
of V = 2af. The scan amplitude was chosen to be between a few beam diameters and an amplitude of one-half the 
diameter of the cloud, about 80 /im. For thermal clouds, we typically scanned the full diameter of the cloud. A 
schematic of the experimental set-up, including the imaging system, is depicted in Fig. |l|. 



Exposing a Bose condensate to the stirring laser beam transferred energy to the atoms, resulting in an increase 
in the thermal component. We measured this component by turning off the scan, allowing the gas to equilibrate for 
100-200 ms and then shutting off the magnetic trap. After ballistic expansion for a fixed time, typically 50 ms, the 
atoms were exposed to a ~ 500 /is pulse of near-resonant light, and the shadow cast by the atoms was imaged onto 
a CCD camera. The two-dimensional transmission profile yielded the column density of the cloud and the kinetic 
energy distribution of the atoms. For mixed clouds this distribution was bimodal, with the condensate localized in the 
center and the more dilute thermal component appearing in the wings. Our technique allowed us to measure energy 
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changes as little as 10 nK, in a regime where the thermal fraction (< 10%) was barely visible in the images. 

We employed this technique in earlier work to analyze the heating induced by stirring the condensate with a blue- 
detuned laser beam at various velocities [0 . There the temperature increase suggested two regimes of dissipation 
separated by a velocity threshold. However, a precise determination of the critical velocity was not possible since 
the sensitivity of the measurements was limited by fluctuations in the background heating processes such as inelastic 
collisions. As a result, the small heating rate due to the stirrer near the critical velocity, of the order of 20 nK/s, 
could not be detected. Moreover, the initial thermal fraction could not be suppressed below 40% due to non-adiabatic 
effects during decompression of the magnetic trap. 

By optimizing the evaporation strategy, we have achieved nearly pure (> 90%) condensates as initial conditions 
for the stirring experiments. While background heating processes could not be fully eliminated, we could account for 
them by subtracting the energy of stirred and unstirred clouds. Thus, even with a nearly invisible thermal component 
we can detect changes in the non-condensed fraction of a few percent. Using this method, we have measured the small 
heating present near the critical velocity that was previously indiscernible. 
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FIG. 2. High sensitivity calorimetry. Shown is an absorption image of an expanded cloud. On the right is a one-dimensional 
vertical slice taken through the center of the image which represents the column density h — — In T, where T is the transmission 
of the probe light. The vertical axis is expanded to see the extremely dilute thermal wings. The solid line is a constrained 
two-dimensional bimodal fit, which yields a thermal fraction of (6 ± 0.6) %. 



Fig. 1^ shows a normalized image of the probe light transmitted by the atomic cloud. The data can be described by 
a bimodal velocity distribution comprised of a Thomas-Fcrmi profile for the condensate and a Gaussian distribution 
for the thermal component M , 
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where p and z are the radial and axial directions, respectively, pc and pT the radial size of condensate (Thomas- Fermi 
radius) and thermal component, respectively, and zc and zt the same for the axial sizes, he and hx are the peak 
column density of the condensate and thermal cloud in the time-of-flight absorption image. The Bose function is 

5n(a:) = Er=i^Vfc"- 

The temperature of the thermal cloud is determined from the width of the Gaussian distribution, 
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where the latter approximation holds for long time of flight ujzt >> 1. The number of thermal atoms is iV^ = 
J hxip, z)dz = Tr[g3{l) / g2ii)]hTPTZT — 2.3 hxpTZT- Due to the high optical density in the center of the images, 
the number of condensate atoms was more reliably extracted from the mean-field energy rather than the integrated 
column density. In the Thomas-Fermi limit this is — ^{15Cj'^h'^\/MNca)'^^^ , where uj = {ujfujzY^^ is the mean 
angular trapping frequency, and a — 2.75 nm is the two-body scattering length for sodium. Equating mean-field 
energy to the measured kinetic energy we get fi — \M{p'^ -\- ^z^)/t^, where t is the time-of-fiight. The total number 
is N = Nt + Nc- 

The transmitted fraction of the probe light was determined by the ratio of images taken with and without the atoms 
present. The two-dimensional data were fitted to the function f{p, z) = Ce~°'"'--'''^\ where a is the known absorption 
cross-section. In general, one can fit all of the unknown parameters Nt/N, Nc /N, N, px, zt, Pc, zc, as well as the 
rotation angle of the image in the p-z plane and a scale factor C caused by fluctuations in the probe light intensity 
between the two images. 
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For very low temperatures the spatial extent Pt,Zt of the thermal component as well as its column density hr 
diminish considerably. It became very difficult to reliably extract thermal fractions below 10% from an unconstrained 
bimodal fit of the form given in Eqn. |[ Therefore, we constrained the thermal fit parameters Nt/N,pt and zt to 
reflect the correlation between the number of thermally excited atoms and the temperature, a relationship which may 
be parameterized in different ways. For example, for an ideal Bosc gas, the thermal fraction and energy per particle 
are 
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where Tc is the transition temperature. Including the effect of the finite mean-field energy fi increases the thermal 
fraction (and therefore the total energy), yielding a correction to the above expressions 
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where 77 = /iT=o/^B7c- We parameterized the functions ^ and r which had been numerically evaluated for specific 
values of 77 in the range 0.3-0.45 ||l5| using mean-field theory. By inverting equation ^ we obtain for the total energy 
of the gas 
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where the function e is determined by the functions ^ and t. The parameters nT,pT, and zt can now be expressed 
in terms of a single parameter, the thermal fraction Nt/N. Without this constraint, the fitting routine tended to 
increase the spatial extent pt, zt of the thermal cloud when the column density ut is low, an unphysical result. 

By imposing this constraint we could extend the bimodal fitting technique to temperatures close to the chemical 
potential p,. Fig. ^ shows a one-dimensional vertical slice taken through the image data, converted into column 
density of the cloud. The thermal wings appear at about 0.3 mm from the center of the image, near the edge of the 
condensate. Although the thermal cloud is extremely dilute, the constrained fitting routine has no difficulty detecting 
the 6% normal fraction of the gas present in the image. With this technique, the energy per particle of the gas can 
be established with an uncertainty of about 10 nK, for temperatures as low as ksT/p — 1.2. To our knowledge, this 
combination of high sensitivity and low temperature thermometry has not been previously reported. 
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FIG. 3. Calorimetry of a condensate. Shown is the energy transfer rate versus velocity of the stirring laser beam measured 
using subtraction of background heating processes (full circles) and older data taken from ref. without background sub- 
traction (crosses). The velocity is in units of the peak sound velocity in the condensate. The error bars reflect shot-to-shot 
variations in the temperature, which is the largest source of error. The inset is a magnification of the region near zero velocity, 
showing sensitivity to ~ 20 nK/s heating rates. 

Fig. H shows the heating of the condensate versus the velocity of the laser beam. The velocity was varied by 
changing the frequency of the laser beam while keeping the scan amplitude fixed at roughly 1 /3 of the axial Thomas- 
Fermi diameter. The heating rate was measured by adjusting the stirring time between 30 ms and 8 s while keeping 
the total energy transfer roughly the same for all velocities used. Earlier data taken by varying the amplitude while 
keeping the frequency fixed at 56 Hz, about 3 times the axial trapping frequency Jill , are shown for comparison. 

The new data show a substantial improvement in signal-to-noise with respect to the older set. This enhanced 
sensitivity allows for a clear resolution of the small heating rate (~ 20 nK/s) visible at low velocities near vjcs — 0.15. 
This is shown in the inset, which is a magnified plot of the data in the low velocity region. These data suggest that 
the critical velocity for excitation occurs at Vc — O.lCs, as confirmed by direct measurements of the drag force from 
the asymmetry of the condensate density profile induced by the scanning laser beam [ p3[ . This technique was shown 
to be consistent with the calorimetric measurements, but had a higher sensitivity at low velocities. 

Above the critical velocity the condensate experiences a drag force F = k(v — Vc), where k is the coefficient of drag 
p6| , p^ . The rate of energy transfer is F ■ v, and so we can find the critical velocity from a parabolic fit above Vc, 

dE 

— = KV[V - Wc) (8) 

This yields a critical velocity of vd Cs — 0.07 ± 0.01 for the new data. This new value is smaller than estimated in our 
earlier paper Both experiments are consistent, since the scatter in the observed heating rates in prevented a 
clear observation of the threshold, and the critical velocity was estimated by linear extrapolation from high heating 
rates to be Vc/cs = 0.25. Other possible extrapolations tend to lower the value of Wc, for instance a function cx v{v — Vc) 
as suggested by equation^, yielded Vc/cg — 0.20 ± 0.07. With our improved calorimetry we have done measurements 
using the fixed frequency /variable amplitude method employed earlier and observed a small heating rate of 20 nK/s 
&t V = 0.15cs that was previously indiscernible. Furthermore, small differences in the laser beam profile may have 
also contributed to the higher value of Vc in . 

Both data sets show good agreement in the overall heating rate. For the new data above v/cg — 0.1 the heating 
rate increases, until about 0.2, where it appears to level off and then to drop to a minimum at around v/cs — 0.25. 
For higher velocities the heating rate increases once again. The suppression of the heating occurring at v/cg — 0.25 
arises from a frequency dependent feature, as discussed in [p^ . At that velocity the laser scans near the axial trapping 
frequency and excites synchronous dipole motion of the condensate. This results in a lower relative velocity between 
condensate and stirrer, thereby reducing the heating. 



IV. ENERGY DISSIPATION IN A NORMAL GAS 



At finite temperature, one would expect heating even below the critical velocity due to physical processes occurring 
within the normal component. For a dilute Bose gas at temperatures below 0.45 Tc, the thermal fraction is less 
than 10%, and plays a very minor role in the heating discussed so far. In addition, the repulsive mean field of the 
condensate lowers the density of the normal component in the center of the cloud even further. This allowed us to 
study the breakdown of superfluidity due to processes initiated primarily within the condensate. 

There is a a major difference in the description of the condensate and the normal component. Interactions within 
the condensate give rise to quantum fluid properties. The macroscopic wavefunction = ^Jne^'^ yields hydrodynamic 
equations of motion for this fluid, expressed in terms of the density n and superfluid velocity v — (?i/m)V(/i. In 
contrast, a dilute normal gas interacts very weakly with itself and with the condensate. For the most part, the mean 
free path for collisions l^i-p > c?, where d is some characteristic length scale, for example, the laser beam diameter or 
the radius of the cloud. In this coUisionless regime, no correlations between particles exist, and a fluid description in 
terms of a local density n and a local velocity v does not hold. However, heating still occurs as long as the thermal 
atoms collide with the moving laser beam. We present measurements of this heating mechanism in a gas above the 
transition temperature, which we can understand using a kinetic model. These measurements allow us to compare 
the relative effectiveness of heating the normal and superfluid components. 
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Let us first consider a one-dimensional model of a hard wall moving at a velocity v through a gas of atoms with 
thermal velocity V ~ ^J^k^TjM (see Fig. ||). We assume v <SiV and elastic collisions between the atoms and the 
wall. 




FIG. 4. One-dimensional model for the heating of a thermal cloud. Atoms from both sides collide at velocity V with a wall 
moving at velocity v. 

An atom moving in the direction of the wall slows down during the collision, while an atom moving against the wall 
increases its velocity. This leads to a net heating rate due to the finite velocity v. The energy transferred to atoms 
on the left and right is 

ASieft = \m{-Y + 2vf - \^MV^ = 2Mv{v - Y) (9) 
Ai^right = \M{y -I- 2vf - ^MV^ = 2Mv{v + V) (10) 
while the corresponding collision rates are 

Fieft = nA{V - v) (11) 
Tright = nA{V + v) (12) 

where n is the atomic density and A the surface area of the moving wall. The rate of energy increase is therefore 

dE 

— = FieftASieft + Fright Alright - SMnAVv^ (13) 

which can be simply interpreted as an energy transfer per collision, cx Mu^, multiplied by the collision rate nAV. The 
model is easily extended to three dimensions and objects of arbitrary shapes. We approximate the stirring laser as an 
infinite cylinder of diameter d moving transversely to its axis of symmetry, along the axial direction of the cylindrically 
symmetric trap. The frequencies of oscillation along the axial and radial direction are ujz and ujr, respectively, and 
the cloud 1/e^ radii are i?^,^ = ^2kBT/Muj'^ ^. Summing over the Boltzmann velocity distribution, and accounting 
for the three-dimensional nature of the collisions, the energy transfer rate per particle is: 

= 2nru'^,Mv^, (14) 



dE 



N 



where the geometrical overlap factor rjA = Ai,/S ^ 2d/TTRr is the ratio between the cross-sectional area of the laser 
beam, Ai,, and that of the thermal cloud, S = ttR^. 
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FIG. 5. Heating of the normal gas. The energy transfer rate per particle versus velocity of the moving laser beam is compared 
to the result of the kinetic model described in the text (solid line). 

Eqn. ^has a simple interpretation: each atom moves back and forth in the axial direction at a rate Vz, and within 
each trapping period a fraction rjA of them hit the stirring beam acquiring an energy oc Mv^ . Fig. ^ shows the energy 
transfer rate to a gas of atoms slightly above the transition temperature, obtained using the calorimetric technique 
described earlier. We varied the velocity by keeping the scan amplitude fixed to approximately the diameter of the 
thermal cloud, and varying the scan frequency. The temperature was obtained from a fit to the wings of the cloud 



[g[ . To compare with the model calculations Eqn. 14 was evaluated for the average of initial and final temperatures 
(the temperature increased by about a factor of 2 during the stirring) using the experimental parameters d = lO/im, 
laser power 18 mW, — 20 Hz. The laser power was chosen to give a barrier height U k, QkgT, where T is the gas 
temperature. The data clearly show the parabolic dependence on the laser beam velocity up to the thermal velocity 
V = jM = 40 mm/s, as predicted from Eqn. [ij, thus verifying the basic tenets of the kinetic model of 

heating. 
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FIG. 6. Dependence of the heating on the power of the stirring laser beam. Shown is the heating of a thermal cloud close 
to the BEG transition temperature versus the laser beam potential Uo in units of the temperature fcsT. Data were taken by 
varying the laser power at a fixed stirring velocity of 7mm/s. The solid line is a fit to the data at low laser power U < lOfcsT. 

The observed absolute heating rate also agreed well with the prediction given by Eqn. |lj to within 20-30%. The 
major uncertainty came from the laser beam size, which may have had a 20-30% variability from run to run. In Fig. 
^we show the heating rate for a cloud near 0.78/iK, plotted versus the ratio of barrier height Uq to temperature ksT, 
where T is the average of initial and final temperatures. In this regime, the heating rate increased in proportion to 
the effective laser beam size, which scales as -y/ln [Uo/kBT). The solid line is a fit to this functional form. At higher 
laser power the data deviate from this fit probably due to additional heating from the non-Gaussian wings of the laser 
beam. 
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FIG. 7. Temperature dependence of the heating. Plotted is the heating coefficient «: described in the text, for two different 
velocities, 9.7 mm/s (x) and 19.4 mm/s (•), versus the average temperature of the gas. The solid line is the value predicted by 
Eqn. The decrease in heating efficiency versus temperature indicates decreasing overlap between the stirrer and the gas. 

In the above model, the temperature of the gas only appears through the factor t^a- Therefore, at higher tempera- 
tures, where the cloud is more spatially extended, one expects the heating rate to decrease as the size Rr increases, or 
= Kw^, where n oc d/^/T. We keep the ratio of barrier height to temperature, and therefore the effective laser 
beam size d, approximately constant. Fig. ^ shows the coefficient of heating k measured for two different stirring 
velocities v below the thermal velocity, along with the model prediction. The data show decreasing efficiency of 
heating at high temperatures, which we attribute to the decrease of geometrical overlap between the cloud and the 
laser beam. 
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FIG. 8. Heating of a thermal cloud close to the BEG transition temperature versus stirring frequency. 

One may ask if the heating measured depends strictly on velocity, or if there is an additional dependence on the scan 
frequency. For low scan frequencies / < 20^^^ we observed no difference in the heating for different frequencies at the 
same laser beam velocity, taking care to avoid the trap resonance at / = i^^. The data presented thus far depended 
only on the laser beam velocity. However, at very high scan frequencies / > 50j^^, the back-and-forth motion of 
the optical potential becomes time- averaged and appears as a quasi-static perturbation to the atoms. Therefore the 
heating rate diminished, as seen in Fig. |^, where the heating rate was measured at a fixed scan amplitude and variable 
frequency. 



V. DISCUSSION 



Since our first report on critical velocities in a condensate, several theoretical papers have added further insight 
^^p3|. The observed critical velocity Vc — O.lCs is lower than the predictions for homogeneous 0,|2^ and inhomo- 



geneous 25 1 2D systems. This discrepancy is most likely due to the actual 3D geometry, where the laser beam pierced 
lower density regions of the condensate. It was noted in Ref. ||2l| that the critical velocity for phonon excitation is 
lowered by the inhomogeneous density distribution. Vortex stretching and half-ring vortices can lower the 3D critical 
velocity below the 2D value . Jackson et al. performed 3D simulations obtaining a critical velocity as low 

as 0.13cs, quite close to our results. 

The relevant critical velocity is most likely related to vortex nucleation |16 r^JT9| , p^ - p6| , which is usually smaller 



than the Landau and Feynman critical velocities ||2^] at which phonons ]2^ , |2l| or vortices |29| , p0t become energetically 
favorable. In the trapped quantum gases, due to the isolation from the environment there is generally no initial 
vorticity present, and vortex nucleation, rather than the energetics of vortex formation, plays the dominant role in 
determining the critical velocity. This is different from the case of liquid ^He, where there is usually vorticity present 
at the surface. As predicted by the Feynman criterion, observed critical velocities depend strongly on the size of the 
flow channel [ p7| which reflects the dependence of the vortex energy on the finite size geometry. 

There are differences between a constantly moving object and our scanning laser beam. The rapid turnaround of 
the laser beam occurring twice per cycle can emit phonons causing additional heating [ p"8|j2^ ]. This heating should 
scale with u^/, where v is the velocity of the scan and / the scan frequency. Our new high sensitivity measurement 
technique could be used to explore this heating mechanism. 

We may compare the superfluid heating with that of the normal gas. Clearly, different mechanisms have been 
outlined for the thermal cloud and for the condensate. For a thermal gas, the heating rate is 



dE 



= Nthtv^ (15) 
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while the condensate heating rate above the critical velocity is 

= Ncncv{v - Vc) ~ NcKcv^ (16) 

c 

where the approximation holds for v 3> Wc- In general we may write the heating rate as A-F, where A is the energy per 
collision with the laser beam, and F ~ nAV is the rate of such collisions. For a thermal gas, V is the thermal velocity 
Vt = y/'^ksT/M. We rewrite this as ^ = (^zRz and use this relation to obtain the velocity V for the condensate as 
Vc = V^Cg, where Rz is the Thomas- Fermi radius. This could also be derived from the pressure nMc^/2 exerted by 
the condensate on the stirrer [^0| . Therefore, in comparing heating in a thermal cloud and a condensate we should 
assume that the condensate atoms strike the stirrer with a velocity Vc oc Cg- The ratio r of the energy transferred 
per collision with the laser beam for normal and condensed gases is then 

At Ntkt ncAcVc ktRt 

r = = • ~ (.1 ' ) 

Ac NcKc tltAtVt kcRc 

where Rt and Rc are the transverse sizes of the normal and condensed components, respectively. We measured 
this ratio as r ~ 2. Thus outside the superfluid regime, the difference in intrinsic heating for normal and superfluid 
components is seen to be quite small. Each collision with the stirrer transfers an energy of about Mv'^ to the cloud, 
both in the thermal cloud and the condensate. Thus the observed difference between the heating coefficients k-t and 
Kc is mainly geometrical due to the different cloud sizes. 

In conclusion, we have studied dissipation induced by a macroscopic stirrer in a Bose gas below and above the BEC 
transition temperature. We present high sensitivity calorimetric measurements showing that the critical velocity is 
lower than suggested by earlier data |1^. Measurements on dilute thermal clouds support a kinetic model of heating. 
This provides a first approximation to dissipation below the BEC transition temperature where one may separately 
consider the contributions from condensed and non-condensed phases. 
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